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Abstract 


The widths for pionic decays of lowest two excited doublets (0+,1*) and 
(1*,2+) of heavy mesons are studied with QCD sum rules in the leading 
order of heavy quark effective theory. The ambiguity due to presence of two 


distinct 1* states are solved. The soft pion approximation is not used. 
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I. INTRODUCTION 


The properties of excited heavy mesons have attracted attention in recent years. The 
mass and decay widths have been studied with potential model combining with data for K 
mesons in [ll], with relativistic Bethe-Salpeter equation in JB] and with QCD sum rules 
in GJ. In work PB] the masses of low lying excited heavy mesons (2*,1*) and (1+, 0*) were 
studied with QCD sum rules in the heavy quark effective theory (HQET) [Å] up to the 
order of O(1/mg), where mg is the heavy quark mass. The present work is devoted to the 
decay widths of these mesons in the same formalism. 

One problem encountered in studying the decay widths of excited heavy mesons with 
QCD sum rules is the following. Except for the lowest states 07, 17, the spectra contains 
a pair of states for any spin-parity j? with close values in their masses but quite different 
in magnitudes of their decay widths. In our case, one of the two 1* states is a narrow 
resonance decaying mainly by emitting a D wave pion, while the other one is a very wide 
resonance decaying by emitting a S wave pion. An interpolating current used for the 
narrow 1* state with a small coupling to the other 1* state may cause sizable error in the 
result of calculation. It is only in the mg — oo limit, there is a conserved quantum number 
je, the angular momentum of the light component, which can be used to differentiate the 
two states. Therefore, HQET has important advantage for this purpose. 

The proper interpolating current J a ae for the states with the quantum number j, P, 


je in HQET was given in [5]. They were proved to satisfy the following conditions 


OA (0) |", P.G = friGyoppo yn , (1) 
: aya; terb . a pull 
i (O|T (75; (2) Jin pr gi (0)) |0) = 605i OPPO ji (1) Sgi PL Oh g 
x f dtåte— vt) pla) (2) 


in the mg > oo limit, where 79% is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, g?? = g*% — v%v% is the transverse metric tensor, 8 


denotes symmetrizing the indices and subtracting the trace terms separately in the sets 


(ai: aj) and (01 --0;), fp;, and IIp,;, are a constant and a function of x respectively 
which depend only on P and jr. Because of equations (1) and (B), the sum rule in HQET 
for decay widths derived from a correlator containing such currents receive no contribution 
from the unwanted states with the same spin-parity as the states under consideration in 
the mg — oo. Starting from the calculations in the leading order, the decay amplitudes 
for finite mg can be calculated unambiguously order by order in the 1/mg expansion in 


HQET. 


II. SUM RULES FOR DECAY AMPLITUDES 


In the present work we shall confine ourselves to the lowest lying excited states in the 
leading order of 1/mg expansion. Denote the doublet (1*,2*) with je = 3/2 by (Bi, B3) 
and the doublet (07, 1+) with jp = 1/2 by (Bj, Bj). From covariance and conservation of 
the angular momentum of the light component in the mg — oo limit, the amplitudes for 
the decay of Bı, Bž to the ground states B, B* in the doublet (07,17) with jp = 1/2 have 


the following forms 


* * v 1 v * 
M(Bi > Ba) = I ein (ata — 394 )9 (Bi, B*) , (3) 
M(B} — Br) = I nur 4 g(B3, B) , (4) 
M(B} > B*r) = I t€apop "na? qug(Bz, B") , (5) 


where Nav, Mu and €, are polarization tensors for states 2*, 1* and 17 respectively. qiu = 
Qu — ty Å = v2, 1 for charged and neutral pion respectively. It can be shown lal 
by combining heavy quark symmetry and chiral symmetry that the constants in (§)-(§) 


satisfy 


9(Bz, B) = g(B3, B*) = E g(Bı, B*) =G. (6) 


For deriving the sum rules for the constants in the decay amplitudes (B), @) and (B) 


we consider the correlators 


= f f atea'y Horr (38 OR ge) 0 
= (apa? ~ sor ?) Go (w,w") , (7) 
fate mar (39000102) 0) = (araf — 39a?) Gazal.) (8) 


[ate eT (77 OIE) 0) = 3 (vgod? +a 8) Gazalo), (9) 


2,+,2 
where k = k' +q, w = 2u-k, w! = w - k' and q? = 0. PCAC is used in the first equality of 
(M. The forms of the right hand side of (f), (H) and (Q) are determined by that a, @ and 
are transverse indices, x — y = vt on the heavy quark propagator and that the pion is in 
the D wave state due to conservation of angular momentum of the light component. The 


interpolationg currents are given in [B] as 


a 1 
ja = 5 ; a a 
Hai [Era (—i) (2; — 3% Pr) q, (10) 


Q1,a 1- (1 ana ama 2 aa 
Tee = I RD (yepp pop — Sop.) a, (11) 


Q 1 L Qa a 1 L 
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where h, is the heavy quark field in HQET and Ytu = Yu — Upp. 


Let us first consider the function Gg, g+ (w, w’) in (M). As a function of two variables, it 


has the following pole terms from double dispersion relation 


DEE ON EE E (13) 
(2A 1 — w’) (2A, 3 — w) 2A 1 — wu! 2A, 3 =W ? 
where fpj, are constants defined in (I), Ap;, = Mpje — mg. As explained in Section [], 
only one state with j? = 1* contributes to as the result of equation (I). This would 
not be true if the last term in (0) is absent. 

Usually, the soft-pion approximation w = w’ is used in QCD sum rules for decay 
amplitudes. This is not a good approximation since the difference w — w’ = Ww -q is 


about 1 GeV in the real decay process as can be seen from the observed spectra of D 


system and the calculation results of [fh]. We instead put 


3 


w-w=%A, 3 -Å 1)=2A. (14) 


2 3 


This is satisfied by the real decay process in the leading order of 1/mg. We then use 
dispersion relation of one variable to Gg, px (4,4), the pole terms of which are 


f_af,,39( BiB") cHe 
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A The interaction of 


The free heavy quark propagator in HQET is fa (x — vr) n 
the heavy quark with the gluon field A in the leading order of 1/mg expansion is ghv - Ah. 
We shall use the Fock-Schwinger gauge x, A(x) = 0 throughout this paper. In this gauge 
the full propagator of the heavy quark (0|T(h(0)A(x))|0) in the leading order of 1/mg is 
identical to the free one, because v = x/7 on the heavy quark line. Therfore, the correlator 


(7) has the form 


får ci" te{P(a(a)|Dea(vr)a(0)10)} = (ata — gi") far f(wg,7), (16) 


where I" is a vertex factor containing y matrices. Since in our approach v -q is fixed to 
be equal to A by (14) and the invaraint function G(B,B*) is considered as a function of 
a single variable w (or w’), the integral in 7 is dominated by small 7, that is, the small 
distance x = v -T for w in the deep euclidian region. Therefore for deriving QCD sum rules 
we can calculate the correlator ({7) by the operator product expansion. This is an expansion 
over the small quantity 7, whch after the integration becomes an assymptotic expansion 
in powers of 1/w. By dimensional reason it is a power series in v - g/w and Agcp/w 
which is meaningful for sufficiently large value of w. After Borel transform the series is 
convergent when v - q/T is sufficiently small, where T is the Borel transform parameter. 
We emphasize that this is so because v -q is fixed in our approach. The perturbation terms 
do not contribute since they do not violate the conservation equation 0" j = 0. For the 
higher terms we expand the quark field q(x) in the matrix element (7(q)|T(q(x)q(0))[0). 
The matrix elements appearing in the expansion have the general form (B7). (BY) and (BY 


in Appendix. The terms in (B7)-(B9) which contribute to the QCD sum rules in the soft 
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pion approximation were used in ÆJ. Here we keep all terms in these formulas in the 
calculations in order to avoid the soft pion approximation. Substituting these equations in 
the correlator (fø) and making Borel transform to supress the contribution of higher states 


we find the sum rule in the leading order of HQET. 


B,B*) så , er 
NE alts (ep cer 4... = 
3 2 2A/T 
59 + ple — da — e2A)Je (17) 


In the above equation we have neglected all radiative corrections and operators with di- 
mension > 6. In the derivation of above equation we need to rotate the integration path 
clock-wise to the imaginary axis of 7 and redefine the integral variable 7' = ir. The result 
is then obtained by using the Borel transform formula BY fl dr er" = (å) H (1/T)” in 
the negative region of w. 

The dots in the left hand side of (17) represents terms from the contributions from the 
radial excited states of the same j?” and jẹ and the continuum. In order to estimate these 
terms we make the duality assumption that above w = w, the dispersion integral of the 
hadron side is equal to that obtained from the QCD calculation. The spectra function 
PQED(x) of the latter is obtained by performing the Borel transform Br to the right 


hand side of (17). The result for pq0p(x) contains two terms proportional to 6(w — 2A) 
œ PQED(W')du' 


; vanishes for we > 2A ~ 1 
We Ww =w 


and 0'(w — 2A) respectively. This implies that 
GeV. Therefore we shall neglect the contributions from the radial excited states and the 
continuum to the sum rule E). 

We can multiply (7) by exp(2A 1 /T) and differentiate in 1/T in order to eliminate c 


and c’. By using (6) we obtain the following sum rule for the constant G 


_ 2A 3 Å JA 
GF fng = HĀ, go + (or — da - eA) (1 + å) a (18) 


Q 


The sum rules for the correlators (H) and (H) can be obtained by similar calculations. 
Equations similar to (I8) are obtained and we have checked that the relations in (  ) are 


satisfied. 


The decay amplitudes of the members of the doublet (j, +, 1/2), Bj and Bj, have the 


following form 


M(B, — Bry =1G', (19) 


M(B, — Btn) = I č - 7G’ 20 
1 


in the mg — co limit. The sum rule in HQET for the coupling constant G’ can be derived 


similarly. From [Å], the interpolating currents for Bj, Bi can be chosen as 
I «= 1 = 
ij = ja = 5a 
Jod = WA hug , Hd = V2 hv. (21) 
As pointed out in [B], for any state (j, P, je) besides the interpolating current with the 
lowest dimension there is another one which has an additional factor —iP, inserting before 
q. For the doublet (j, +, 1/2), they are 


f 1 => 1 1 po 
f — ; i CEN 5n, a i 
Joni ar; h(-i) Pig , Fed 2 hyn i) Peg - (22) 


As discussed in PB], there is some motivation for using the latter currents for the corre- 
sponding constant fia j2 Survives in the non-relativistic limit of the Bethe-Salpeter wave 
function, while the constant f;.1/2 for the current (21) vanishes in this limit, since Bj and 
B; are P wave states in the non-relativistic quark model. 

After Borel transformation and eliminating the single pole terms by differentiation, the 


sum rule derived by using the interpolating current is 


fattig = 2[A; 1(—goA? + 3h2A + 3h1) — (—3b1 + 31A + 3d;A 


2 
2A i ogi 
+ dog? + 2dA? — AL Få AT, (23) 


where A = M 12 —Å- 1. If the interpolating current (2J) is used instead of (22), the 


sum rule turns out to be 


= A 18—7y ÅL yy og T 
OE (Ayo -E a a +3/ 


where 


a= 2(fi+ fed), B= 4(gA + pA? + hi), 


y = A(ay + ag Å + 2c,A + eA? + er A3) (25) 


and A= Axe = Aaz 


III. DETERMINATION OF THE PARAMETERS 


In order to obtain the decay constants G and G' from (8), (23) and (24) we need to use 
the mass parameters A’s and the coupling constants f’s of the corresponding interpolating 
currents as input. Å 1/2 and f 1/2 can be obtained from the results in [SB] as Å 1/2 =0.5 
GeV and f 1/2 0.24 GeV? at the order a, = 0. Notice that the coupling constant -1/2 
defined in the present work is a factor 1/v2 smaller than that defined in (Al. Ais 123 Aga /2 
and Nå j2 are given in BJ. fi,3/2, fi,1/2 and f} 1/2 can be determined from the formulas 
(34), (27) and (28) of reference [5] derived from sum rules for two point correlators. The 


results are 


Å, 3/2 = 0.95 GeV f+3/2 = 0.19 +0.03 Gev*? , 
Mj. =0.9 GeV Fap = 0.284 0.03 GeV*? , 
Ayia = 1.05 GeV fia =0.3240.03 GeV?” . (26) 


(40)-(45) still contains several parameters. må is known to be 0.8 GeV’. mî was 
calculated in P] with QCD sum rules with the result m? = 0.2 GeV’. We have obtained 
QCD sum rules for ej and e3 by considering the two point functions 


b 


i | e@ate(0lT{a(v)ins sale) 400) 5 tn D)*insa(0)|0) = (n-gPAW@)bn, (27) 
i f etatz (Ol {ala yrs ale) (0) 5 (n- DY? hsg(0)0) = (na), (29 


where n? = 0. The resulting QCD sum rules are 


2e = 1 PO /M? 1 2 112 (qa) 
= = oe ds + — (aG TOs ; 29 
z, T Tr) ET ET: 29) 


7 


3 80 m2 
2e2 fr må = aa | (m4 — sje %/M" ds F Fay os Yo (: + za) 
20 må \ (Gay 
+ 2710 2 + A) M2 ` (30) 


In obtaining (BO) we have eliminated the a; pole term by differentiation with respect to 
the Borel parameter 1/M?. ma œ 1.3 GeV is the mass of aq. In the numerical calculations 


for e and e we use the following standard values of the vacuum condensates: 


(Gq) = —(0.24 GeV)? , 


(asGG) 9 = 0.038 GeV* (31) 


as well as a, = 0.4. From (29) and (BO) we obtain the following numerical values for e1 


and ez, 


ej = —0.014 + 0.002 GeV’, 


eg = 0.015 + 0.002 GeV. (32) 


IV. NUMERICAL RESULTS AND DISCUSSION 


We now turn to the numerical evaluation of the sum rules for G and G’. By using (26) 
and (Bll) together with the values of parameters determined above we obtain G and G' as 
functions of T from the sum rules (18), (23) and (24). The results are plotted as curves in 
fig. 1. The lower limit of T is determined by the requirement that in the sum rules before 
differentiating with respect to 1/T (such as (L), the terms of higher order in 1/T in the 
operator expansion are reasonably smaller than the leading term, say < 1/3 of the latter. 
This leads to T > 1.2 GeV and T > 3 GeV for the sum rule (18) and (23) respectively. 

Usually the value of T is constrained from the upper side by the requirement that the 
contribution from higher states is sufficiently suppressed. The continuum model discussed 
below equation ([7) allows T to be any large value. However, in order that the results 


depend less on the continuum model it is good if T is not much larger than the 2(A“) — A), 


8 


where mg + A® is the mass of the radial excited state. Quark model calculations 
indicate that AY — A is roughly 0.7 GeV. Therefore, this requirement can be satisfied for 
the sum rule for G. We find within the stability window subject to the constrain disussed 


above 
G=5.7+0.6+0.7 GeV’, 13 GeV <T <1.5 GeV. (33) 


From the view point stated above, the sum rule (23) for G’ is not very satisfactory since 
the lower limit of T is too high. If we somewhat arbitrarily take the width of the window 
to be 0.3 GeV, we find from (23) 


G'=10+0.1+£0.5, 3 GeV <T <3.3 GeV. (34) 


In (BJ) and (34) the first error refers to the variation with T in this region, while the second 
error takes into account the uncertainty in the values of the coupling constant f’s. The 
errors quoted above do not include those from higher order terms in OPE and the 1/mg 
expansion and the values of parameters characterizing the condensations. 

However,the (83) and (B4) have another source of uncertainty. These sum rules depend 
on the parameter ej which is obtained by the sum rule (29). It is well known [P] that the 
sum rules derived from the correlator of two pseudo-scalar currents give too small value for 
the physical quantities. In the case of two pseudo-scalar current without derivative, the 
phenomenogical side is about a factor 2 larger than the theoretical side. Just for a rough 
estimation, let us multiplying the ej value by a factor 2. We find that the T value at which 
the last term is 1/3 of the first term in the sum rule before differention with respect to 
1/T move to lower values 0.6 GeV and 1.9 GeV for the sum rule (17) and (RJ) respectively. 


That is, the convergence becomes faster. The values of G and G' change to 
G=36410+40.7 GeV”, 1.0 GeV<T<1.5 GeV, 


G'=12+401+405, 18 GeV <T<21 GeV. (35) 


On the other hand the sum rule (24) is not good in stability. Especially, G’ changes 


sign near T œ~ 1.2 GeV. Therefore, we do not use it in the numerical calculation of G’. 
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The decay widths in the leading order are 


P(B: > Ba) = "lal, 

I(B; > Br) = 0'li , 

(Bj — 8"n) = la, 
P(B; > Br) = 29 


3 y 
(Bi > B'n) = 20 ld, (36) 


where sum over charged and neutral pion final states has been included in the above 
formulas. 

In the following we apply the leading order formulas obtained above to the excited states 
of charmed mesons. The value g for the first three processes in (BG) are calculated from the 
experimental mass values of the relevent particles and that for the last two processes 
are calculated from the difference between the theoretical mass values of the two doublets 
in the leading order obtained in [Å] and [B]. For the excited states of D we obtain the 
results listed in table fl, where the numbers in the parentheses are results for e} = —0.028. 
The experimental data [1] and the results of BB] obtained with relativistic B-S equation are 
also included. For the D series the O(1/mg) corrections are not negligible. Yet, assuming 
the total width of D> is dominated by one pion decay, the calculation result for D3 is in 
reasonable agreement with experimental data. The calculation result for I (D; — D*7) in 
the leading order is significantly smaller than the experimental data for [(D, — all). This 
can partly be attributed to mixing of the two 1* states in the order of O(1/mg) [4]. This 
effect can be calculated within our formalism by introducing the O(1/mo) terms in the 
Largrangian of HQET. In reference [fl], it was estimated that the [(D; > Dp + D*p) ~ 6 
MeV, ['(D3 — Dp+ D*p) ~ 3 MeV. Therefore, the effect of mixing of two 1 states in the 
order O(1/mg) needs not to be very large. The calculated results for the decay widths of 


Dy and D1 in the leading order is too small. This may due to the slow convergence of Eq. 
(B3). 
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Appendix: The OPE of matrix element (7'(q)|T(q(x)q(0))|0) 


(DIECO) = ilfrs + fadder (2) (37) 
(x"(q)| Dual (0)45(0)10) = {(91 + gafas + hitus + haier Jas (F) > (88) 
(aa Dy Dea (0)R(0)IO) = ik [guns + aaf) — iopvbi + ee Yabo 
+ (Qu + AYn)er + (Ww = WIn)cr 


+ (Gur + Wor) C2 + UGOrw — Qaru) do 


Ti 
++ edaabhsda (3) : (39) 
1 T 
f= 3 00 i fo = -4 ; (40) 
1 å Le 

gı ap, anno > 2 fe » hj=0, h= 127 (44)o ; (41) 

må 5 1 fr 
a= Gr (ao » Q= — 3g fr > a= går » C= 247, (Aa)o > (42) 


å _ Ti 5 
where (7'(a|a 5 D°y159|0) = —ifr mi dp - 


i i _ Ti 
er(n a)? = ;(m(glaz (n: D)o), for n*=0, (43) 
1 - = Ti 
ea(n- 4) = z(o la z hln D)o), for n° =0, (44) 
a e 1 
bi = = 5 by = —3c 5 di = 0 5 dy = os + 12f, (dq)o ‘ (45) 


Here fr = 93 MeV. 
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Figure Captions 


Fig. 1. Dependence of G and G' on the Borel parameter T. The three 


curves correspond, from top to bottom, to the sum rules (1g), (23) and 


BA, respectively. 
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TABLES 
TABLE I. The decay widths of D1, D3, D) and Dj. The numbers in the parenthesises are 


results obtained by using e; = —0.028 instead of —0.014. 


Transition Calculated Width(MeV) Observed Width(MeV) Result from 
[(D, > D*r) 7.0(2.8) 8.0 
T(D; — all) 18.9+48 

oe 7 

I (Di > D*r) 6.9(2.8) 6.8 
(Dx — all) 23+% 

T(D) > Dr) 

ID > D*7) 45(64) 
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